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O ■ ABSTRACT 

Q ■ The Generalized Chaplygin Gas (GCG) with the equation of state p = was 

^ . recently proposed as a candidate for dark energy in the Universe. In this pa- 

. per we confront the GCG with SNIa data using available samples. Specifically 

>• . we have tested the GCG cosmology in three different classes of models with (1) 

>< : ^m = 0.3, Qch = 0.7; (2) il^ = 0.05, fich = 0.95 and (3) ^rn = 0, Qch = 1, 

as well as a model without prior assumptions on Qrn- The best fitted models 
are obtained by minimalizing the function. We supplement our analysis with 
confidence intervals in the [Aq, a) plane by marginalizing the probability density 
functions over remaining parameters assuming uniform priors. We have also de- 
rived one-dimensional probability distribution functions for Qch obtained from 
joint marginalization over a and Aq. The maximum value of such PDF provides 
the most probable value of Qch within the full class of GCG models. The general 
conclusion is that SNIa data give support to the Chaplygin gas (with a = 1). 
However noticeable preference of Aq values close to 1 means that the a depen- 
dence becomes insignificant. It is reflected on one dimensional PDFs for a which 
turned out to be fiat meaning that the power of present supernovae data to dis- 
criminate between various GCG models (differing by a) is weak. Extending our 
analysis by relaxing the prior assumption of the flatness of the Universe leads to 
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the result that even though the best fitted values of are formally non-zero, 
still they are close to the flat case. Our results show clearly that in GCG cos- 
mology distant (i.e. z > 1) supernovae should be brighter than in ACDM model. 
Therefore one can expect that future supernova experiments (e.g., SNAP) hav- 
ing access to higher redshifts will eventually resolve the issue whether the dark 
energy content of the Universe could be described as a the Chaplygin gas. More- 
over, it would be possible to differentiate between models with various value of a 
parameter and/or discriminated between GCG, Gardassian and AGDM models. 
This discriminative power of the forthcoming mission has been demonstrated on 
simulated SNAP data. 

Subject headings: cosmology: theory — distance scale — supernovae: Generalized 
Ghaplygin Gas 



1. Introduction 

For a couple of years two independent observational programs — the high redshift 
supernovae surveys (Perlmuttcr ct al. 1999) and CMBR small scale anisotropy measurements 
(dc Bcrnardis ct al. 2000, Benoit ct al. 2003, Hinshaw et al. 2003) have brought a new 
picture of the Universe in the large. While interpreted ithin the FRW models results of 
these programs suggest that our Universe is flat (as inferred from the location of acoustic 
peaks in GMBR power spectrum) and presently accelerates its expansion (as inferred from 
the SNIa Hubble diagram) . Combined with the independent knowledge about the amount of 
baryons and GDM estimated to be Vim = 0.3 (Turner 2002) it follows that about Qx — 0.7 
fraction of critical density per = %Jq should be contained in a mysterious component 
called "dark energy". The most obvious candidate for this smooth component permeating 
the Universe is the cosmological constant A representing the energy of the vacuum. Well 
known fine tuning problems led many people to seek beyond the A framework, and the 
concept of the quintessence had been conceived. Usually the quintessence is described in a 
phenomenological manner, as a scalar field with an appropriate potential (Ratra & Peebles 
1988, Galdwell, Dave & Steinhardt 1995, Prieman, Stebbins & Waga 1995). It turns out, 
however, that quintessence program also suffers from its own fine tuning problems (Kolda & 
Lyth 1999). 

In 1904 Russian physicist Ghaplygin introduced the exotic equation of state p — 
to discribe an adiabatic aerodynamic process (Ghaplygin 2004). The attractiveness of this 
equation of state in the context dark energy models comes mainly from the fact that it gives 
the unification of both dark energy (postulated in cosmology to explain current aceleration 
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of the Universe) and clustered dark matter which is postulated in astrophysics to explain 
the fiat rotation curves of spiral galaxies. It is interesting that the Chaplygin gas can be 
derived from the quintessence Lagrangian for the scalar field with some potential and also 
from the Born-lnfeld form of the Lagrangian (Kamenshchik;, Moschella & Pasquier 2001). 
The Chaplygin equation of state has some interesting connections with string theory and it 
admits the interpretation in the framework of brane cosmologies (Jackiw 2000). 

Recently the so called Chaplygin gas (Kamenshchik, Moschella & Pasquier 2001, Fabris, 
Gongalves & de Souza 2002, Szydlowski & Czaja 2004) — a hypothetical component with the 
equation of state p = — was proposed as a challenge to the above mentioned candidates 
for dark energy. This, also purely phenomenological, entity has interesting connections with 
string theory (Ogawa 2000). Currently its generalizations admitting the equation of state 
p — —-^ where < a < 1 have been proposed (Bento, Bertolami & Sen 2002, Carturan & 
FineUi 2002a). 

In this paper we confront the Generahzed Chaplygin Gas with the SNIa data. At this 
point our choice of Generahzed Chaplygin Gas cosmologies deserves a sort of justification. 
There are two approaches in the literature. First one is phenomenological, namely having no 
preferred theory of dark energy responsible for acceleration of the Universe one characterizes 
it as a cosmic fluid with an equation of state px = wpx where w > —1 (see e.g. (Chiba 
1998, Turner & White 1997) and an immense literature that appeared thereafter). Because, 
as already mentioned above, a strain of ideas about dark energy is associated with an evolving 
scalar there arc good reasons to expect that cosmic equation of state could be time dependent 
i.e. w = w{t) = w{z) (e.g. Weller & Albrecht 2001, Maor et al. 2001 and many others 
thereafter). This approach seems attractive from the perspective of analyzing observational 
data such like supernovae surveys and indeed this approach was taken while first analyzing 
the data (Perlmutter et al. 1999, Knop et al. 2003 or Riess et al. 2004). However even though 
such analysis places constrains on any potential theory that might explain the dark energy 
phenomenon, ultimately one always ends up at testing a specific theory. Along this line there 
appeared attempts to reconstruct the scalar potential, if the scalar field was responsible for 
dark energy (e.g. Alam et al. 2003 and references therein). Our approach goes along this 
philosophy but instead is devoted to the Generalized Chaplygin Gas which is being recently 
considered as candidate to unified dark matter-energy component (i.e. responsible for both 
clustering and accelerated expansion (Makler, de Oliveira & Waga 2003). 

The cosmological models with the Generalized Chaplygin Gas have also many special 
features which make them atractive. In standard cosmological model one can clearly distin- 
guish the epochs of radiation domination followed by (ordinary) matter domination (with 
decelerated expansion). As mentioned above supernova data suggest that the epoch of de- 
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celerated expansion ended and switched to accelerated epoch — dominated by dark energy. 
The Generahzed Chaplygin Gas models describe smoothly the transition from the deceler- 
ated to accelerated epochs. They represent the simplest deformation of concordance ACDM 
(Gorini et al. 2004). And moreover, they propose a new unified macroscopic (phenomenolog- 
ical) description of both dark energy and dark matter. This places them in a distinguished 
position from the point of view of Occam's razor principle. It should be also noted that the 
Genarahzed Chaplygin Gas model allowed us to explain presently observed acceleration of 
the Universe without the cosmological constant and/or modification of Einstein's equations. 

If one takes seriously given dark energy scenario (necessary to explain cosmic acceler- 
ation) one should also consider the behaviour of perturbations in such a universe. In the 
framework of quintessence models with the barotropic equation of state (i.e. p = wp and 
w — const) one faces the problem of instabilities in short scales. This appears because the 
speed of sound squared (equal here to w) is negative (and constant). Calculation of the sound 
speed in Generalized Chaplygin Gas model (see below) reveal its non-barotropic nature. The 
perturbations in GCG models are stable in short scales even in an accelerating phase (Car- 
turan & Finelli 2002a). Moreover, they behave like dust perturbations when Chaplygin Gas 
is in dust regime. 

Another motivation for studying Generalized Chaplygin Gas models goes from theo- 
retical physics — specifically from attempts to describe the dark energy in terms of the 
Lagrangian for a tachyonic field (Garousi 2000, Sen 2002). Of course it would be nice to 
have a description of dark energy in terms of the non quintessence Lagrangian as it describes 
the nature of dark energy while the cosmological constant is only phenomelogical and ef- 
fective description. One should also note that the Generalized Chaplygin Gas equation of 
state arises in modern physics in the context of brane models (Bordemann & Hoppe 1993, 
Kamenshchik, Moschella & Pasquier 2001, Randall & Sundrum 1999) where the Generalized 
Chaplygin Gas manifests itself as an effect of immersion of our Universe in multidimensional 
bulk space. 

Generalized Chaplygin gas models have been intensively studied in the literature and in 
particular they have been tested against supernovae data (Makler, de Ohveira & Waga 2003, 
Avclino et al. 2003, Collistete et al. 2003), Icnsing statistics (Dev, Alcaniz & Jain 2003), 
CMBR measurements (Bento, Bertolami & Sen 2003a, 20003b, Carturan k Finelh 2003b, 
Amendola et al. 2003), age-redshift relation (Alcaniz, Jain & Dev 2003), x-ray luminosities of 
galaxy clusters (Cunha, Lima & Alcaniz 2003) or from the large scale structure considerations 
(Bean & Dore 2003, Multamaki, Manera & Gaztanaga 2003, Bilic et al. 2003). Perspectives 
to distinguish between Generalized Chaplygin Gas, brane-world scenarios and quintessence 
in forthcoming gravity wave experiments has been discussed in (Biesiada 2003). Although 
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the results are in general mutually consistent there was no strong convergence to unique 
values of o; parameters characterizing Chaplygin gas equation of state. 

Makler, de Ohveira & Waga (2003) have considered the FRW model filled completely 
with Generahzed Chaplygin Gas and concluded that whole class of such models is consistent 
with current SNIa data although the value oi a — 0.4 is favoured. This result has been 
confirmed by our analysis (class (3) models). However, when the existing knowledge about 
baryonic matter content of the Universe was incorporated into the study our results were 
different from Makler, de Oliveira & Waga (2003) who found that a = 0.15 was preferred 
(assuming flrn = 0.04 which is very close to our assumption for class (2) models). 

As noticed by Bean & Dore (2003) Generalized Chaplygin Gas models have an inherent 
degeneracy with cosmological constant models as far as background evolution is concerned, 
and therefore they have a good fit with SNIa data. These degeneracies disappear at the 
level of evolution of perturbations and hence confrontation with CMBR spectrum would be 
decisive. Using available data on the position of CMBR peaks measured by BOOMERANG 
(de Bernardis et al. 2000) and Archeops (Benoit et al. 2003, Hinshaw et al. 2003, Bento, 
Bertolami & Sen (2002) obtained the following constraints: 0.81 < Aq < 0.85 and 0.2 < a < 
0.6 at 68% CL in the model representative of our class (2) (i.e. with = 0.05 assumed). 
Another estimation of the parameter a was done by Amendola et al. (2003) with WMAP 
Data. The obtained the < a < 0.2 at 95% confidence level. 

Using the angular size statistics for extragalactic sources combined with SNIa data it 
was found in (Alcaniz & Lima 2003) that in the the 0,m — 0.3 and Qch — 0.7 scenario 
best fitted values of model parameters are Ao = 0.83 and a = 1. respectively. Recent 
paper by Bertolami et al. (2004) in which Generalized Chaplygin Gas models have been 
analyzed against Tonry et al. (2003) supernovae data relaxing the prior assumption on 
fiatness suggests, surprisingly as the authors admit, the preference of a > 1. 
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Table 1. Results of statistical analysis of Generalized Chaplygin Gas model (with 
marginalization over A4) performed on analyzed samples of SNIA (A, C, K6, K3, TBI, 

TBll, Silver, Gold) as a minimum best-fit (denoted BF) and with the maximum 
likelihood method (denoted L). First two rows for each sample refer to no prior on Q^- 
The same analysis was repeated with fixed priors = 0.0, Qrn = 0.05 and = 0.3. 



sample 




^Ch 


Ao 


a 


M 




method 


A 


0.00 


1.00 


0.77 


1.00 


-3.39 


95.4 


BF 




0.17 


0.83 


0.83 


0.00 


-3.36 


— 


L 




0.00 


1.00 


0.77 


1.00 


-3.39 


95.4 


BF 




0.00 


1.00 


0.73 


1.00 


-3.38 


— 


L 




0.05 


0.95 


0.80 


1.00 


-3.39 


95.4 


BF 




0.05 


0.95 


0.76 


1.00 


-3.38 


— 


L 




0.30 


0.70 


0.96 


1.00 


-3.39 


95.8 


BF 




0.30 


0.70 


0.96 


0.00 


-3.38 


— 


L 


C 


0.00 


1.00 


0.80 


1.00 


-3.44 


52.9 


BF 




0.15 


0.85 


0.86 


0.00 


-3.41 


— 


L 




0.00 


1.00 


0.80 


1.00 


-3.44 


52.9 


BF 




0.00 


1.00 


0.76 


0.49 


-3.43 


— 


L 




0.05 


0.95 


0.83 


1.00 


-3.44 


53.0 


BF 




0.05 


0.95 


0.79 


0.11 


-3.43 


— 


L 




0.30 


0.70 


0.99 


1.00 


-3.42 


53.3 


BF 




0.30 


0.70 


0.99 


0.00 


-3.39 


— 


L 


K6 


0.00 


1.00 


0.81 


1.00 


-3.52 


55.3 


BF 




0.10 


0.90 


0.88 


0.00 


-3.51 


— 


L 




0.00 


1.00 


0.81 


1.00 


-3.52 


55.3 


BF 




0.00 


1.00 


0.78 


0.71 


-3.52 


— 


L 




0.05 


0.95 


0.84 


1.00 


-3.52 


55.4 


BF 




0.05 


0.95 


0.81 


0.06 


-3.52 


— 


L 




0.30 


0.70 


1.00 


1.00 


-3.51 


55.9 


BF 




0.30 


0.70 


1.00 


0.00 


-3.49 


— 


L 


K3 


0.00 


1.00 


0.85 


1.00 


-3.48 


60.4 


BF 




0.11 


0.89 


0.88 


0.00 


-3.45 


— 


L 




0.00 


1.00 


0.85 


1.00 


-3.48 


60.4 


BF 




0.00 


1.00 


0.80 


0.30 


-3.47 




L 




0.05 


0.95 


0.87 


1.00 


-3.47 


60.4 


BF 




0.05 


0.95 


0.84 


0.00 


-3.47 




L 




0.30 


0.70 


1.00 


1.00 


-3.44 


61.4 


BF 




0.30 


0.70 


1.00 


0.00 


-3.42 




L 


TBI 


0.00 


1.00 


0.79 


1.00 


15.895 


273.9 


BF 




0.00 


1.00 


0.81 


1.00 


15.905 




L 




0.00 


1.00 


0.79 


1.00 


15.895 


273.8 


BF 




0.00 


1.00 


0.75 


1.00 


15.905 




L 




0.05 


0.95 


0.82 


1.00 


15.895 


274.0 


BF 




0.05 


0.95 


0.78 


1.00 


15.915 




L 




0.30 


0.70 


0.97 


1.00 


15.915 


275.8 


BF 




0.30 


0.70 


0.96 


0.00 


15.915 




L 


TBII 


0.00 


1.00 


0.78 


1.00 


15.915 


186.5 


BF 




0.00 


1.00 


0.81 


1.00 


15.925 




L 




0.00 


1.00 


0.78 


1.00 


15.915 


186.5 


BF 




0.00 


1.00 


0.75 


1.00 


15.915 




L 




0.05 


0.95 


0.81 


1.00 


15.915 


186.6 


BF 




0.05 


0.95 


0.78 


1.00 


15.925 




L 
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2. Cosmological model 

Einstein equations for the Friedman-Robertson- Walker model with hydrodynamical 
energy-momentum tensor T^^, = (p + p)Ujj,Ui, — pg^^ read: 

8nGp k 



a J 3 a'^{t) 

d(t) A-kG , „ , 

^ = -— (p + 3rt (2) 

Let us assume that matter content of the Universe consists of pressure-less gas with en- 
ergy density p„j representing baryonic plus cold dark matter (CDM) and the Generalized 
Chaplygin Gas with the equation of state 

representing the dark energy responsible for the acceleration of the Universe. If one further 
makes an assumption that these two components do not interact, then the energy conserva- 
tion equation 

p + ?,H{p + p) = (4) 

where H = a/ais the Hubble function, can be integrated separately for matter and Chaplygin 
gas leading to well known result pm = PmfiO'^^ and (see also Bento, Bertolami & Sen 2002 
or Carturan & Finelli 2002) 

/ B \^ 

The physical interpretation of, so far arbitrary, constants A and B is the following. 
Adopting usual convention that current value of the scale factor gq is equal to 1, one can 
see that pch,o = {A + B)^ represents the current energy density of the Chaplygin gas. 
Calculating the adiabatic speed of sound squared for the Chaplygin gas 

9 dpch OiA aA 

c: = 



dpcH P'^'^ A + 



it is easy to confirm that the current value of cj. is equal to c^ q = -^f^- Hence the constants A 
and B can be expressed as combinations of quantities having well defined physical meaning. 

Our further task will be to confront the Chaplygin gas model with SNIa data and for 
this purpose we have to calculate the luminosity distance in our model 
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Table 1 — Continued 



sample 




^Ch 


Ao 


a 


M 




method 




0.30 


0.70 


0.97 


1.00 


15.925 


188.4 


BF 




0.30 


0.70 


0.96 


0.00 


15.935 




L 


Silver 


U.UU 




n Qo 
U.oz 


1 nn 
l.UU 


io.y40 




J3r 




0.00 


1.00 


0.84 


1.00 


15.945 




L 




0.00 


1.00 


0.82 


1.00 


15.945 


229.4 


BF 




0.00 


1.00 


0.79 


1.00 


15.955 




L 




0.05 


0.95 


0.85 


1.00 


15.945 


229.6 


BF 




0.05 


0.95 


0.81 


1.00 


15.955 




L 




0.30 


0.70 


0.99 


1.00 


15.965 


232.3 


BF 




0.30 


0.70 


0.99 


0.00 


15.965 




L 


Gold 


0.00 


1.00 


0.81 


1.00 


15.945 


173.7 


BF 




0.00 


1.00 


0.83 


1.00 


15.955 




L 




0.00 


1.00 


0.81 


1.00 


15.945 


173.7 


BF 




0.00 


1.00 


0.77 


1.00 


15.955 




L 




0.05 


0.95 


0.84 


1.00 


15.945 


173.8 


BF 




0.05 


0.95 


0.80 


1.00 


15.955 




L 




0.30 


0.70 


0.99 


1.00 


15.965 


175.6 


BF 




0.30 


0.70 


0.99 


0.00 


15.965 




L 



Table 2: Generalized Chaplygin Gas model parameter values obtained from the marginal 
probability density functions calculated on Perlmutter, Knop, Tonry/Barris and Riess sam- 
ples with prior relaxed. 



sample 








^Ch 






a 


A 


0. 


1 7+0.08 
' ' -0.17 


0, 


oq+0.17 

'°'-'-0.08 


0. 


oq+0.14 
'°'-'-0.09 


_0.0+0.67 


C 


0. 


1 p:+0.08 
'-'-3-0.15 


0. 


or+0.15 
'°'-'-0.08 


0. 


of^+0.13 
■^O-O.IO 


0.0+0.66 


K6 


0. 


-IQ+O.ll 

■J^'J-o.io 


0. 


qO+0.10 


0. 


OO+0.12 
■OO-0.08 


_0.0+0.66 


K3 


0. 


-1 -1 +0.07 
■^^-0.11 


0. 


oq+o.ii 

■°^-0.07 


0. 


oo+O.ll 
■OO-0.05 


0.0+0.66 


TBI 


0. 


00+0.21 


1. 


.00_o.21 


0. 


81 +0-12 
■"-'--0.07 


1-0-0.60 


TBII 


0. 


00+0.21 


1. 


.00_o.21 


0, 


81 +0-12 
'°-'--0.07 


1-0-0.62 


Silver 


0. 


00+0.18 


1, 


.00_o.l8 


0, 


8/1+0-09 


1-0-0.59 


Gold 


0. 


00+0.20 


1. 


.00_o.20 


0. 


oo+O.ll 
■°'^-0.07 


1-0-0.64 
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Table 3: Generalized Chaplygin Gas model parameter values obtained from the marginal 
probability density functions calculated on Perlmutter, Knop, Tonry/Barris and Riess sam- 
ples. The analysis was done with fixed O^n = 0.0, = 0.05 and ilm = 0.3. 



sample 




^Ch 


Ao 


a 


A 


0.00 


1.00 


0.73t°:?« 


1-0-0.63 




0.05 


0.95 




1-0-0.66 




0.30 


0.70 


0.96+°:°^ 


00+0.65 


C 


0.00 


1.00 


76+0-08 


40+0-36 
'-'•+^-0.35 




0.05 


0.95 


0.79+H? 


41+0-27 




0.30 


0.70 


0.99t°:?} 


0,0+0.64 


K6 


0.00 


1.00 


0.78^H^ 


71+0-29 
' J.-0.40 




0.05 


0.95 


n oi +0.08 

U.5i_o.o9 


'"'•'J '-'-0.06 




0.30 


0.70 


1.00_o.io 


00+0.64 


K3 


0.00 


1.00 




0-3ol°:i 




0.05 


0.95 


o.84l:H? 


00+0.67 




0.30 


0.70 


1.00_o.06 


00+0.63 


TBI 


0.00 


1.00 


' '-'-0.05 


1-0-0.54 




0.05 


0.95 


n 7O+0.04 


1-0-0.55 




0.30 


0.70 


qfi+0-04 


00+0.67 


TBll 


0.00 


1.00 




1-0-0.54 




0.05 


0.95 




1-0-0.54 




0.30 


0.70 


0.961H^ 


00+0.67 


Silver 


0.00 


1.00 




1-0-0.52 




0.05 


0.95 


A 0-1 +0.04 
U-8-'--0.04 


1-0-0.54 




0.30 


0.70 


n qq+0.01 


00+0.64 


Gold 


0.00 


1.00 


0.77tH^ 


1-0-0.58 




0.05 


0.95 


o.solH^ 


1-0-0.59 




0.30 


0.70 


0.99lg;°l 


00+0.64 
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Table 4: Results of statistical analysis of Generalized Chaplygin Gas models with flat prior 
relaxed and with marginalization over Ai performed on Knop Samples K3 as a minimum 
best-fit (denoted BF) and with the maximum likelihood method (denoted L). First two rows 
refer to no prior on Q„i. The same analysis was repeated with fixed = 0.0, Q^n = 0.05 



and i}„i 


= 0.3. 
















sample 






^Ch 


^0 


a 


M 




method 


K3 


-0.19 


0.00 


1.19 


0.82 


1.00 


-3.48 


60.3 


BF 




-0.60 


0.00 


1.26 


0.89 


0.00 


-3.46 




L 




-0.25 


0.00 


1.25 


0.82 


1.00 


-3.49 


60.3 


BF 




0.10 


0.00 


0.90 


0.76 


0.00 


-3.46 





L 




-0.28 


0.05 


1.23 


0.84 


1.00 


-3.49 


60.3 


BF 




0.05 


0.05 


0.90 


0.78 


0.00 


-3.47 




L 




-0.48 


0.30 


1.18 


0.93 


0.97 


-3.49 


60.3 


BF 




-0.35 


0.30 


1.05 


0.88 


0.00 


-3.47 




L 


Gold 


-0.12 


0.00 


1.12 


0.80 


0.99 


15.945 


173.4 


BF 




-0.32 


0.00 


1.06 


0.82 


0.00 


15.945 




L 




-0.13 


0.00 


1.13 
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Table 5: Results of statistical analysis of Generalized Chaplygin Gas models with fiat prior 
relaxed and with marginalization over Ai performed on Knop Samples K3. Model parameter 
values are obtained from the marginal probability density functions. First row refer to no 
prior on Q^. The same analysis was repeated with fixed fi^ = 0.0, flm = 0.05 and flm = 0.3. 
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Fig. 1. — Residuals (in mag) between the Einstein-de Sitter model (zero line), the ACDM 
model (upper curve) and the best-fitted Generalized Chaplygin Gas model with = 
0.3, flch = 0.7 (middle curve), sample K3. 
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Fig. 2. — Levels of constant on the plane {AQ,a) for Generalized Cliaplygin Gas model 
with flrn = 0.3, flch = 0.7, sample K3, marginalized over 7W. The figure shows preferred 
values of Aq and a. 
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Fig. 3. — Confidence levels on the plane {Ao, a) for Generalized Chaplygin Gas model with 
Qjn. = 0.3. Qch = 0.7, sample K3, marginalized over A4. The figure shows the ellipses of 
preferred values of Aq and a. 
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Fig. 4. — Residuals (in mag) between the Einstein-de Sitter model (zero line), the flat 
ACDM model (upper curve) and the best-fitted Generalized Chaplygin Gas model with 
= 0.05, flch = 0.95 (middle curve), sample K3. 
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Fig. 5. — Levels of constant on the plane {AQ,a) for Generalized Chaplygin Gas model 
with flrn = 0.05, flch = 0.95, sample K3, marginalized over A4. The figure shows preferred 
values of and a. 
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Fig. 6. — Confidence levels on the plane {Aq, a) for Generalized Chaplygin Gas model with 
Qm = 0.05, Qch = 0.95, sample K3, marginalized over A4. The figure shows the ellipses of 
preferred values of Aq and a. 
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Fig. 7. — Residuals (in mag) between the Einstein-de Sitter model (zero line), the flat 
ACDM model (upper curve) and the best-fitted Generalized Chaplygin Gas model with 
= 0, flch — 1 (middle curve), sample K3. 
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Fig. 8. — Levels of constant on the plane {AQ,a) for Generalized Chaplygin Gas model 
with Qm = 0, flch = 1, sample K3, marginalized over A4. The figure shows preferred values 
of Aq and a. 
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Fig. 9. — Confidence levels on the plane {Aq, a) for Generalized Chaplygin Gas model with 
Qm = 0, Qch = 1, sample K3, marginalized over A4. The figure shows the ellipses of preferred 
values of Aq and a. 




Fig. 10. — The density distribution (one dimensional PDF) for Qch obtained from sample K3 
by marginalization over remaining parameters of the model. We obtain the limit D,ch > 0.70 
at the confidence level 95.4%. 
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where — —wi ^^^^ 

J-{x) = sinh{x) 
J-'{x) = X 
J-{x) = sin{x) 



for k<0 (7) 
for k = (8) 
for k>0 (9) 



The Friedman equation (1) can be rearranged to the form giving exphcitly the Hubble 
function H{z) = a/a 



H{zf = Hi 



\lm{l + Zf + ^Ch {Ao + (1 - A)(l + Zf^^^'^^) + Qfe(l + Zf'\ (10) 



where the quantities Vti, i = m,Ch,k represent fractions of critical density currently con- 
tained in energy densities of respective components and Qm + ^ch + — 1- For the 
transparency of further formulae we have also denoted Aq = A/{A + B). 

Finally the luminosity distance reads: 

c 1 



dL{z) = {l + z) 

-^0 Vl^fel 



dz' 



+ zf + Qch + (1 - + + 0^(1 + z)\ 



(11) 



The formula (11) is the most general one in the framework of Friedman- Robertson- 
Walker cosmology with Generalized Chaplygin Gas. Please note, that this model propose a 
unified macroscopic (phenomenological) description of both dark energy and dark matter. 

Further in this paper we will mostly use its version restricted to flat model k — Q (the 
exception will be when we relax flat prior) since the evidence for this case is very strong in 
the light of current CMBR data. Therefore while talking about model testing we actually 
mean the estimation of a and Aq parameters for the best fitted flat FRW cosmological model 
filled with Generalized Chaplygin Gas. 

To proceeded with fitting the SNIa data we need the magnitude-redshift relation 

m{z, M, l^m, ^ch, Aq.q) = M + h logio Dl{z, VL^, VLch] Aq, a) (12) 

where: 

Dl{z, Qrn, ^Ch', Aq, a) = H^diiz, Ho, Urn, ^Ch', Aq, a) 
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is the luminosity distance with Hq factored out so that marginahzation over the intercept 

A4 = M-51ogio^o + 25 (13) 

leads actually to joint marginahzation over Hq and M (M being the absolute magnitude of 
SNla). 

Then we can obtain the best fitted model minimalizing the function 



where the sum is over the SNIa sample and ai denote the (full) statistical error of magnitude 
determination. This is illustrated by figures (Fig.l - Fig.9) of residuals (with respect to 
Einstein-de Sitter model) and levels in the (Aq, a) plane. One of the advantages of 
residual plots is that the intercept of the m — z curve gets cancelled. The assumption that 
the intercept is the same for different cosmological models is legitimate since Ai is actually 
determined from the low-redshift part of the Hubble diagram which should be linear in 
all realistic cosmologies. From the full Perlmutter's sample (see below) we have obtained 
A4 = —3.39 which is in a very good agreement with the values reported in the literature 
(Perlmutter et al. 1999, Riess et al. 1999). For other samples see discussion below. 

The best-fit values alone are not relevant if not supplemented with the confidence levels 
for the parameters. Therefore, we performed the estimation of model parameters using 
the minimization procedure, based on the likelihood function. We assumed that supernovae 
measurements came with uncorrelated Gaussian errors and in this case the likelihood function 
jC could be determined from chi-square statistic C oc exp (— x^/2) (Perlmutter et al. 1999, 
Riess et al. 1999). 

Therefore we supplement our analysis with confidence intervals in the (^40, a) plane by 
calculating the marginal probability density functions 

V{Ao,a)(x J exp{-x\nrn,nch,Ao,a,M)/2)dM 

with Qrn: ^Ch fixcd (Qm — 0.0, 0.05, 0.3 ) and 

V{Ao,a)(x / ex.p{-x^{nm,^Ch,Ao,a,M)/2)dnm 



with M. fixed {A4 — —3.39) respectively (proportionafity sign means equal up to the nor- 
malization constant) . In order to complete the picture we have also derived one-dimensional 
probability distribution functions for ftch obtained from joint marginahzation over a and 
Aq. The maximum value of such PDF informs us about the most probable value of flch 
(supported by supernovae data) within the full class of Generalized Chaplygin Gas models. 
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Fig. 11. — The density distribution (one dimensional PDF) for Aq obtained from sample K3 
by marginalization over remaining parameters of the model. We obtain the limit Aq > 0.74 
at the confidence level 95.4%. 




Fig. 12. — The density distribution (one dimensional PDF) for a obtained from sample K3 
by marginalization over remaining parameters of the model. We obtain the limit a < 0.94 
at the confidence level 95.4% . 
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Fig. 13. — Residuals (in mag) between the Einstein-de Sitter model (zero line), the flat 
ACDM model ( upper curve) and the best-fitted Generalized Chaplygin Gas model (without 
any prior assumptions on Qrn) (middle curve), sample K3. 
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Fig. 14. — Levels of constant on the plane (Aq, a) for Generalized Chaplygin Gas model, 
sample K3, marginalized over A4, and Q^. The figure shows preferred values of and a. 
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Fig. 15. — Confidence levels on the plane {Aq, a) for Generalized Chaplygin Gas model, 
sample K3, marginalized over A4, and flm- The figure shows the ellipses of preferred values 
of Aq and a. 
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Fig. 16. — Levels of constant on the plane (Q^, a) for Generalized Chaplygin Gas model, 
sample K3, marginalized over A4, and Aq. The figure shows preferred values of Q^, and a. 
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Fig. 17. — Confidence levels on the plane {Qm, a) for Generalized Chaplygin Gas model, 
sample K3, marginalized over M., and ^40. The figure shows the ellipses of preferred values 
of f^TO^and a. 
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Fig. 18. — Residuals (in mag) between the Einstein-de Sitter model (zero line), flat ACDM 
model (two upper curves: for SNIA with z < 1 — curve located higher and for all supernovae 
la belonging to the sample — curve located lower) and the best-fitted Generalized Chaplygin 
Gas model (without any prior assumptions on Q^) (middle curve), GOLD sample. 
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Fig. 19. — Levels of constant on the plane {Aq, a) for Generalized Chaplygin Gas model, 
Gold sample, marginalized over M., and Q^. The figure shows preferred values of Aq and a. 
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Fig. 20. — Confidence levels on the plane {Aq, a) for Generalized Chaplygin Gas model, Gold 
sample, marginalized over J\4, and flm- The figure shows the ellipses of preferred values of 
Aq and a. 
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3. Fits to ^0 and a parameters 

3.1. Samples used 

Supernovae surveys (published data) have already five years long history. Beginning 
with first published samples other data sets have been produced either by correcting original 
samples for systematics or by supplementing them with new supernovae (or both). It is not 
our intention here to suggest a distinguished role to any one of these data sets. Therefore, 
in our analysis we decided to use a collection of samples from all existing supernovae data. 

Latest data was compiled by Riess et al. (2004). However, for comparision and ilus- 
tration we analyse three other main samples of supernovae. Such investigation seems to be 
useful because it is pointed out in the literature that result of the analysis with different 
SNla sample give often different results (see for example Godlowski Szydlowski & Krawiec 
2003, Choudhury & Padmanabhan 2004). 

Samples from the original Perlmutter et al. (1999) data chosen for the analysis comprise 
the full sample reported by Perlmutter (sample A) and a sub-sample after excluding two 
outliers differing the most from the average lightcurve and two outliers claimed to be likely 
reddened (sample C). Although the outliers often suggest statistical inhomogeneity of the 
data (and some hints suggesting the necessity of removing them from sample A exist) there is 
always a danger that removal of outliers is to some extent subjective. Therefore we retained 
the full sample A in our analysis. The Perlmutter data were gathered some years ago, hence 
one should also refer to more recent supernovae data as well. 

Recently Knop et al. (2003) have reexamined the Permutter's data with host-galaxy ex- 
tinction correctly assessed. Prom the Perlmutter's sample they chose only these supernovae 
which were spectroscopically safely identified as type la and had reasonable color measure- 
ments. They also included eleven new high redshift supernovae and a well known sample 
with low redshift supernovae. In Knop et al. (2003) a few subsamples have been distin- 
guished. We considered two of them. The first is a subset of 58 supernovae with corrected 
extinction (Knop subsample 6; hereafter K6) and the second is that of 54 low extinction 
supernovae (Knop subsample 3; hereafter K3). Samples C and K3 are similarly constructed 
as containing only low extinction supernovae. The advantage of the Knop sample is that 
Knop's discussion of extinction correction was very careful and as a result his sample has 
extinction correctly applied. 

Another sample was recently presented by Tonry et al. (2003) who collected a large 
number of supernovae data published by different authors and added eight new high redshift 
SN la. This sample of 230 SN la was re-calibrated with a consistent zero point. Wherever 
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possible the extinction estimates and distance fitting were recalculated. Unfortunately, one 
was not able to do so for the full sample (for details see Tabic 8 in Tonry et al. (2003)). This 
sample was further improved by Barris et al. (2003) who added 23 high redshift supernovae 
including 15 at 2; > 0.7 thus doubling the published record of objects at these redshifts. 
Despite of the above mentioned problems, the analysis of our model using this sample of 
supernovae could be interesting. Hence for comparison, we decided to repeat our analysis 
with the Tonry /Barris sample. 

We decided to analyze two Tonry/Bariss subsamples. First, we considered the full 
Tonry/Barris sample of 253 SNc la (hereafter sample TBI). The sample contains 218 SNc la 
with low extinction. Because Tonry's sample has a lot of outliers especially at low redshifts, 
we decided to analyze the sample where all low redshift {z < 0.01) supernovae were excluded. 
In the sample of 193 supernovae all SN la with low redshift and high extinction were removed 
(hereafter sample TBII). 

Tonry et al. (2003) and Barris et al. (2003) presented the data of redshifts and luminos- 
ity distances for their supernovae sample. Therefore, Eqs. (12) and (13) should be modified 
appropriately (WiUiams et al. 2003): 

m - M = 5 logio(^L)Tonry " 5 log^Q 65 + 25 (14) 

and 

M = -5 logio //o + 25. (15) 

For the Hubble constant Hq = 65 km s~^ Mpc~^ one gets M = 15.935. 

Recently Ricss ct al. (2004) significantly improved the former Ricss sample. They 
discovered 16 new type la Supernovae. It should be noted that 6 of these objects have 
z > 1.25 (out of total number of 7 object with so high redshifts). Moreover, they compiled 
a set of previously observed SNIa relying on large, published samples, whenever possible, 
to reduce systematic errors from differences in calibrations. With this enriched sample it 
became possible to test our prediction that distant supernovae in GCG cosmology should be 
brighter than in ACDM model (see discussion below). This is the reason why we repeated 
our analysis with new Riess sample. 

The full Riess sample contains 186 SNIa ("Silver" sample). On the base of quality of 
the spectroscopic and photometric record for individual Supernovae, they also selected more 
restricted "Gold" sample of 157 Supernovae. We have separately analyzed ACDM model for 
supernovae with z < 1 and for all SNIa belonging to the Gold sample. 
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3.2. Cosmological models tested 

On these samples we have tested Generahzed Chaplygin gas cosmology in three different 
classes of models with (1) = 0.3, Qch = 0.7; (2) Qm = 0.05, Qch = 0.95 and (3) flm = 0, 
flch = 1- We started with a fixed value of = —3.39 modifying this assumption accordingly 
while analyzing different samples. 

The first class was chosen as representative of the standard knowledge of Qm (baryonic 
plus dark matter in galactic halos (Peebles & Ratra 2003)) with Chaplygin gas responsible 
for the missing part of closure density (the dark energy). 

In the second class we have incorporated (at the level of Qm) the prior knowledge about 
baryonic content of the Universe (as inferred from the BBN considerations). Hence this class 
is representative of the models in which Chaplygin gas is allowed to clump and is responsible 
both for dark matter in halos as well as its diffuse part (dark energy). 

The third class is a kind of toy model - the FRW Universe filled completely with 
Chaplygin gas. We have considered it mainly in order to see how sensitive the SNIa test is 
with respect to parameters identifying the cosmological model. 

Finally, we analyzed the data without any prior assumption about 

3.3. Results 

In the first class of models best fit (with fixed value oi A4 — —3.39) from the sample A 
is {a = 1, Ao = 0.96) at the = 95.8. Sample C gives the best fit of {a = 0.95, Ao = 0.95) 
at the = 53.6. 

In the second class Sample A gives the best fit of {a = 1,^0 = 0.80) at the x^ = 95.4 
whereas the sample C gives the best fit {a — 0.51, Aq — 0.73) at the x^ = 53.7. 

Finally, in the third class the sample A again gives the best fit of (a = 1, = 0.77) at 
the x^ = 95.4 whereas the sample C gives the best fit {a — 0.42, = 0.69) at the x^ = 53.7. 

It should be noted, however that the fitting procedure for sample C prefers A4 = —3.44 
instead of Ai = —3.39 as for sample A. If one takes this vahic the results for sample C 
will change respectively and then for the first class Aq — 1 (at — 53.5) what means 
(see equation (10)) that a can be arbitrary and the problem is effectively equivalent to the 
model with cosmological constant. Analogously, for the second class = 0.83 a = 1 (at 
X^ = 52.9), while for third class Aq = 0.80 a = 1 (at = 52.9). This indicates clearly that 
model parameters, especially a strongly depend on the choice of M.. 
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Separately we analyzed the data without any prior assumption about Qjn- For the 
sample A we obtain as a best fit (minimizing x^) = 0., {a = I.Aq = 0.77) at the 
= 95.4. For sample C assuming A4 = —3.39 we obtain = 0.27, {a = 1,Aq = 0.93) 
with = 53.6 while for M = -3.44 the best fit gives = 0., {a = l,^o = 0.80) with 
X' = 52.9. 

Using minimization procedure, based on likelihood function, joint marginalization over 
Aq and a gives the following results. For the sample A we obtain Qch = 0.82 (hence 
flm = 0.18), with the limit Qch > 0.76 at the confidence level 68.3% and flch > 0.69 at the 
confidence level 95.4%. For the sample C we obtain (for A4 = —3.39) Qch = 0.76 (hence 
n^n = 0.24), with the limit Qch(^{0.69, 0.94) on the confidence level 68.3% and Qch > 0.62 at 
the confidence level 95.4%. For M = -3.44 we obtain: Qch = 0.84 (hence fi^ = 0.16), with 
the hmit Qch^{0.79, 0.98) on the confidence level 68.3% and flch > 0.69 at the confidence 
level 95.4%. 

One could see that results are different for different values of the intercept A4 in each 
sample. Therefore we additionally analyzed our samples marginalized over M.. The results 
are displayed in Table 1. First rows for each sample correspond to no prior on assumed. 
Two fitting procedures were used: fitting (denoted as BF) and maximum 1 ikelihood 
method (denoted L). 

Table 2 contains the results of joint marginalization over Q^, Aq and a. For sample 
A we obtain Qch — 0.83 (hence = 0.17), while for sample C we obtain Qch — 0.85 
(hence Qm — 0.15). For and a we obtain values {a — 0,Ao — 0.83) for sample A, while 
{a — 0,Ao — 0.86) for sample C. With the marginalization procedure we can also obtain 
one dimensional probability distribution function (PDF) for Aq and a for particular class 
models with fixed fl^ (Table 3). 

From the above mentioned analysis we concluded that and Aq parameters derived 
from samples A and C are similar. For Qm fixed Aq increases with increasing Q„i- The 
estimates of a parameter are different for each of two above mentioned samples, but unfor- 
tunately errors are big. Discrepancy between the best fitting procedure and minimization 
procedure (based on likelihood function) increases with the number of parameters fitted. 
The minimization procedure seems to be more appropriate in the context of our problem. 

One can see from the Table 1 that using the Knop's samples had not influenced conclu- 
sions in a significant way. However, the errors of parameter estimation decreased noticeably 
(see Tables 2 and 3). Minimization procedure prefers (especialy for sample K3) a close to 
zero. The exception is the model with = where a — 0.3 and a — 0.71 are obtained for 
the samples K3 and K6 respectively. 
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The above mentioned results for the Knop sample K3 are illustrated on figures Fig. 1 - 
17. In Fig.l we present residual plots of redshift-magnitudc relations between the Einstein-de 
Sitter model (represented by zero line) the best-fitted Generalized Chaplygin Gas model with 
Q„i = 0.3, ^Ich = 0.7) (middle curve) and the flat ACDM model with A = 0.75 and = 0.25 
(Knop et al. 2003)) - upper curve. One can observe that systematic deviation between 
ACDM model and Generalized Ghaplygin Gas model gets larger at higher redshifts. The 
Generahzed Chaplygin Gas model predict that high redshift supernovae should be brighter 
than that predicted with ACDM model. 

Levels of constant on the {Aq, a) plane for Generalized Chaplygin Gas model with 
flrn = 0.3, Qch = 0.7), marginalized over Ai arc presented in Figure 2. The figure shows 
preferred values of Aq and a. Figure 3 displays confidence levels on the {Aq, a) plane for 
Generahzed Chaplygin Gas model with = 0.3, Qch — 0.7) , marginahzed over M.. This 
figure shows the ellipses of preferred values of {Aq and a) . Similar results for the models with 
Qrn — 0.05, Qch — 0.95 and = 0, Qch — 1 are presented in the figures 4-9 respectively. 

Separately we repeated our analysis without prior assumptions on Q^. The density 
distribution (one dimensional PDF) for model parameters obtained by marginalization over 
remaining parameters of the model are presented in Figures 10-12. 

Residuals (in mag) between the Einstein-de Sitter model (zero line), flat ACDM model 
- upper curve - and the best-fitted Generalized Chaplygin Gas model (without prior as- 
sumptions on Q^) - middle curve - are presented on Figure 13, while levels of constant 
and confidence levels on the {Aq, a) plane (marginalized over A4) are presented on Figs. 14 
and 15. 

One should notice that as a best fit we obtain = 0, Qch = 1, = 0.85, a — 1 i.e. 
results are the same as for a toy model with Chaplygin gas only {^ch — 1)- Formally, we 
could have analyzed models with a > 1. However, due to large error in estimation of the a 
parameter, it does not seem reasonable to analyze such a possibility with current supernovae 
data. 

Levels of constant the (Q^, oi) plane for Generahzed Chaplygin Gas model, 

marginalized over Ai and Aq arc presented on Figure 16, while confidence levels on the 
{Qm,Oi) plane (marginalized over A4, and Aq) are presented on Figure 17. This figure shows 
that all three model with Qch = 1, ^Ch = 0.95, and flch = 0.7, are statistically admissible 
by current Supernovae data. 

The results of similar analysis obtained with Tonry/Barris sample arc similar to those 
obtained with previous samples. For example TBll sample gives the best fit: flm — 0, flch — 
1, Aq = 0.78, a = 1 i.e. nearly the same as in the case of K3 sample. 
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Joint marginalization over parameters gives the following results: Qch — 1-00 (hence 
O,^ — 0.0), with the limit Qch > 0.79 at the confidence level of 68.3% and flch > 0.67 
at the confidence level of 95.4%. {a = 1.0, = 0.81) with the limit a E (0.40,1.) and 
Ao e (0.74, 0.93) at the confidence level of 68.3% and a e (0.06, 1.) and Aq e (0.70, 1.00) at 
the confidence level of 95.4%. 

However with minimization procedure we find important difference bectwecn results 
obtained with Tonry/Barris sample and those obtained with Perlmutter and Knop sam- 
ples. Minimization procedure (only except the model with fixed = 0.3) performed on 
Tonry/Barris data gives o = l. It is significantly different from the result obtained for C 
Perlmutter 's and Knop's samples where minimization procedure preffered small values of a 
parameter. Also the Tonry /Barris sample preffered value of D,m — while Perlmutter and 
Knop samples sugested is close to zero, what indicates that barionic component is small 
and in agreement with BEN. 

The new Riess sample leads to the results which are similar to these obtained with 
Tonry/Bariss sample. However the errors in estimation of the parameters are lower. For the 
Gold sample, joint marginalization over parameters gives the following results: flch — 1-00 
(hence Qm — 0.0), with the limit flch > 0.80 at the confidence level of 68.3% and Q,ch > 0.69 
at the confidence level of 95.4%. (a = 1.0, = 0.83) with the limit a G (0.36,1.) and 
Ao e (0.76, 0.94) at the confidence level of 68.3% and a e (0.05, 1.) and Aq e (0.72, 1.00) at 
the confidence level of 95.4%. 

As one can see from the Figure 18 the differences between the results obtained in both 
cases arc small (however the result obtained with the full sample leads to the prediction of 
brighter distant supernovae than in the case with z < 1 SNla.) Residuals (in mag) between 
the Einstein-de Sitter model (zero line), fiat ACDM model (two upper curves: for SNIA 
with z < 1 - higher curve - and for all supernovae la belonging to the sample - lower curve) 
and the best-fitted Generalized Chaplygin Gas model (without prior assumptions on Q^) 
(middle curve) are presented on Figure 18. Figure 18 shows that most distant supernovae 
are actually brighter than predicted in ACDM model. This is in agreement with prediction 
of the Generalized Chaplygin Gas cosmology. 

The levels of constant confidence levels on the {Aq, a) plane (marginalized over 

M.) are presented on Figs. 19 and 20. Fig.20 shows that confidence levels on the {Aq, a) 
(for Gold Sample) plane are comparable at the the 95.4% confidence level with the results 
obtained on the Knop's sample. However the preferred values of a are different. 
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3.4. Flat prior relaxed 

We extended our analysis by adding a curvature term to the original GCG model. Then 
in equation (11) we must take into account the Qk term. For statistical analysis we restricted 
the values of the fl^n parameter to the interval [0, 1], flch to the interval [0,2] and flk "was 
obtained from the constraint fl^ + flch + ^k — 1- However, the cases Qjt < — 1 were excluded 
from the analysis. The results are presented in Table IV and V. 

In the model without prior assumptions on we obtain with Knop's sample flk = 
—0.19 as a best fit, while with maximum hkelihood method prefers — —0.60. However, 
the model with priors on or il,ch the maximum likelihood method prefers the universe 
much "closer" to the flat one. Specifically, for the "toy" model with Chaplygin Gas only one 
gets Q,k — 0.10 and Qk — 0.05 for the model with baryonic content only i.e. — 0.05. One 
should emphasize that even though we allowed flk ^ ^ the preferred model of the universe 
is nearly a fiat one, which is in agreement with CMBR data. It is an advantage of our GCG 
model as compared with ACDM model where in Perlmutter et al. (1999) high negative value 
of Qk was obtained as a best fit, although zero value of was statistically admissible. In 
order to find the curvature of the Universe they additionally used the data from CMBR and 
extragalactic astronomy. 

Density distribution functions (one dimensional PDF) for model parameters obtained 
by marginalization over remaining parameters of the model are presented in Figures 21- 
25. For Qfc we obtain the hmit G (—0.98, —0.22) at the confidence level 68.3% and 
Jljfc e (—1,0.23) at the confidence level 95.4%. For a and Aq parameters we obtain the 
following results: a = 0. and Aq = 0.89 with the limit a e (0, 0.64) and Ao e (0.82, 1) at the 
confidence level 68.3% and a e (0.,0.95) and Aq e (0.73, 1.) at the confidence level 95.4%. 
For the density parameter Qch we obtain the limit flch £ (0.87, 1.51) at the confidence level 
68.3% and flch G (0.61, 1.79) at the confidence level 95.4%. For the density parameter flm 
we obtain the limit Qm < 0.29 at the confidence level 68.3% and flm < 0.53 at the confidence 
level 95.4%. 

Our main result is that the preference of the nearly fiat universe is confirmed with the 
new Riess sample. In the model without a prior assumption on Qrn we obtain flk = —0.12 as 
a best fit, with the Gold sample while maximum likelihood method prefers = —0.32 i.e 
the Gold sample gives even "more fiat" universe than Knop's sample. The models with priors 
on give also very similar results when we analyse Knop and Riess samples. One can see 
that estimation of other models parameters give similar result for both samples only with 
exception for parameter a. Specifically, for the "toy" and "baryonic" models the maximum 
likelihood method prefers the universe with non zero parameter a like for "flat universe". 
One can see that, with flat prior relaxed, when we analyse the Gold sample, the errors in 
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estimation of the model parameter significantly decrase with comparision to the case of the 
Knop sample. 

4. Generalized Chaplygin Gas model in perspective of SNAP data 

In the near future the SNAP mission is expected to observe about 2000 SN la supernovae 
each year, over a period of three years ^. Therefore it could be possible to discriminate 
between various cosmological models since errors in the estimation of model parameters 
would decrease significantly. We tested how a large number of new data would influence the 
errors in estimation of model parameters. We assumed that the Universe is flat and tested 
three classes of cosmological models. In the first ACDM model we assumed that = 0.25, 
and = 0.75 and Ai = —3.39 (Knop et al. 2003). Second class was representative of the 
so called Cardassian models (Freese & Lewis 2002) with parameters flm = 0.42, Qcard = 0.52 
and n — —0.77 as obtained in (Godlowski, Szydlowski & Krawiec 2004). Let us note, that 
technically i.e. at the level of tests hke the Hubble diagram Cardassian models are equivalent 
to quintessence models. The difference is in the underlying philosophy: quintessence assumes 
exotic dark energy component with hydrodynamical equation of state in ordinary FRW model 
while the Cardassian Universe assumes modification of the Friedman equation (which can 
be either due to exotic matter component or due to modification of gravity law). The last 
model is Generalized Chaplygin Gas Model with parameters obtained in the present paper 
as best fits for the K3 sample {flm = 0, Aq = 0.85, a = 1.). These values are in agreement 
with results of the analysis performed on Tonry/Barris and Riess samples. Alternatively, we 
also test the Generahzed Chaplygin Gas Model with small value of a parameter suggested 
by analysis of the Perlmmutter and Knop samples {Qm — 0, Aq — 0.76, a — 0.40). For three 
above mentioned models we generated a sample of 1915 supernovae (Samples Xl,X2,X3a,X3b 
respectively) in the redshift range ze[0.01, 1.7] distributed according to predicted SNAP data 
(see Tab I of Alam et al. (2003)). We assumed Gaussian distribution of uncertainties in 
the measurement of m and z. The errors in redshifts z are of order Icr = 0.002 while 
uncertainty in the measurement of magnitude m is assumed la = 0.15. The systematic 
uncertainty is agys = 0.02 mag at z = 1.5 (Alam et al. 2003). Hence one can assume that 
o'sysiz) = (0.02/1.5)2; in first approximation. For such generated sample we repeated our 
analysis. The result of our analysis is presented on figures Fig.26-29. On these figures we 
present confidence levels on the plane {Aq, a) for sample of simulated SNAP data. The 
figures show the eUipses of preferred values of Aq and a. It is easy to see that with the 
forthcoming SNAP data it will be possible to discriminate between predictions of ACDM 



http://www-supernova.lbl.gov, http://snfactory.lbl.gov 
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and GCG models. With Cardasian Model the situation is not so clear, however (see Fig 
27 and 28). Note that if a ~ 0.4 as suggested by analysis of Perlmutter sample C, (see 
also: Makler, de Oliveira & Waga 2003, Avelino et al. 2003, Fabris, Gonngalves & de Souza 
2002, Collistete et al. 2003) then it will be possible to discriminate between model with 
Chaplygin gas and Cardassian model (see Fig 27 and 29). Moreover, it is clear that with the 
future SNAP data it would be possible to differentiate between models with various value of 
a parameter. This is especially valuable since all analyses performed so far have had weak 
sensitivity with respect to a. 



5. Conclusions 

It is apparent that Generalized Chaplygin Gas models have brighter supcrnovac at 
redshifts z > 1. Indeed one can see on respective figures (Figs 1, 4, 7, 13) that systematic 
deviation from the baseline Einstein de Sitter model gets larger at higher redshifts. This 
prediction seems to be independent of analysed sample. 

Wc obtained that the estimated value of Aq is close to 0.8 in all considered models 
with exepction the model class (1) {VLm = 0.3) when Aq > 0.95. Extending our analysis by 
relaxing the flat prior lead to the result that even though the best fitted values of 0,^ are 
formally non-zero, yet they are close to the fiat case. It should be viewed as an advantage of 
the GCG model since in similar analysis of ACDM model in Perlmutter et al. (1999) high 
negative value of Vlk were found to be best fitted to the data and independent inspiration 
from CMBR and extragalactic astronomy has been invoked to fix the curvature problem. 
Another advantage of GCG model is that in a natural way we obtained the conclusion that 
matter (baryonic) component should be small what is in agreement with prediction from 
BBN (Big Bang Nucleosynthesis). Both estimations of Aq, and VLrn are independent of 
the sample used in our analysis. 

Our results suggest that SNIa data support the Chaplygin gas (i.e. a = 1) scenario when 
the best fitting procedure is used. The minimization procedure performed on Tonry/Barris 
and Riess data gives also a = 1 (only except the model with fixed Vim = 0.3). However, the 
maximum likelihood fitting with Knop et al.'s sample prefers, quite unexpectedly, a small 
value of a or even a = 0, i.e. the ACDM scenario. Please note that small value of a is 
in agreement with the results obtained from CMBR (de Bernardis et al. 2000 Benoit et 
al. 2003, Hinshaw et al. 2003, Bento, Bertolami & Sen 2002, Amendola et al. 2003) and 
with the recent analysis of Zhu (2004) who, using combined data of X-ray gas mass fraction 
of the galaxy cluster, FR IIB radiogalaxies and combined sample Perlmutter et al. (1998) 
and Riess et al. (1998, 2001), sugested that a could be even less than 0. The results are 
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dependent both on the sample chosen and on the prior knowledge of A4 in which the Hubble 
constant and intrinsic luminosity of SNIa are entangled. Moreover the observed preference 
of Aq values close to 1 means that the a dependence becomes insignificant (see equation 
(10)). It is reflected on one dimensional PDFs for a which turned out to be flat meaning 
that the power of the present supernovae data to discriminate between various Generalized 
Chaplygin Gas models (differing by a) is weak. 

However, we argue that with future SNAP data it would be possible to differentiate 
between models with various value of a parameter. Residual plots indicate the differences 
between ACDM and Generalized Chaplygin Gas cosmologies at high redshifts. Therefore one 
can expect that future supernova experiments (e.g. SNAP) having access to higher redshifts 
will eventually resolve the issue whether the dark energy content of the Universe could be 
described as a Chaplygin gas. The discriminative power of forthcoming SNAP data has been 
illustrated on respective figures (Fig.26-29) obtained from the analysis on simulated SNAP 
data. 
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Fig. 24. — The density distribution (one dimensional PDF) for Qk obtained from sam- 
ple K3 by marginalization over remaining parameters of the model. We obtain the limit 
Qjfce(-1,0.23) at the confidence level 95.4%. (Non-flat GCG model.) 
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Fig. 25. — The density distribution (one dimensional PDF) for Q^n obtained from sample K3 
by marginalization over remaining parameters of the model. We obtain the limit 0^6(0, 0.53) 
at the confidence level 95.4%. (Non-flat GCG model.) 




Fig. 26. — Confidence levels on the plane {Aq, a) for sample XI (ACDM model) of simulated 
SNAP data, marginalized over Q^. The figure shows the ellipses of preferred values of Aq 
and a. 
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Fig. 27. — Confidence levels on the plane {Aq, a) for sample X2 (Cardassian model) of 
simulated SNAP data, marginalized over fl^. The figure shows the ellipses of preferred 
values of and a. 
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Fig. 28. — Confidence levels on the plane {Aq, a) for sample X3a (Generalized Chaplygin 
Gas model) of simulated SNAP data {Q^ = 0, Aq = 0.85, a — 1.), marginalized over Q^. 
The figure shows the ellipses of preferred values of Aq and a. 
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Fig. 29. — Confidence levels on the plane {Aq, a) for sample X3b (Generalized Chaplygin 
Gas model) of simulated SNAP data {Q^n = 0, Aq = 0.76, a ~ 0.49), marginalized over Q^. 
The figure shows the ellipses of preferred values of and a. 
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